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We study the propagation of pions at finite temperature and finite chemical potential in the 
framework of the linear sigma model with 2 quark flavors and N c colors. The velocity of massless 
pions in general differs from that of light. One-loop calculations show that in the chiral symmetry 
broken phase pions, under certain conditions, propagate faster than light. 
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I. INTRODUCTION 

The linear sigma model, originally proposed as a model 
for strong nuclear interactions pj, today serves as an ef- 
fective model for the low-energy (low-temperature) phase 
of quantum chromodynamics. The model exhibits spon- 
taneous breaking of chiral symmetry and restoration at 
finite temperature. In the chiral symmetry broken phase 
at zero temperature pions, being massless, propagate 
with the velocity of light. At finite temperature one 
would expect chiral pions to propagate slower than light, 
owing to medium effects. Indeed, this expectation has 
been confirmed by one-loop calculations in the linear 
sigma model with no fermions U- Lil- LJ- Lj| ■ 

In this paper we discuss the propagation of pions at 
nonzero temperature and nonzero finite baryon density 
in a model with two quark flavors. It turns out that 
pions in the presence of fermions become superluminal 
in a certain range of temperature and baryon chemical 
potential. 

A superluminal propagation has recently been studied 
in the context of the Casimir effect in quantum electro- 
dynamics. Scharnhorst has demonstrated |fj that, when 
vacuum fluctuations of the electromagnetic field obey pe- 
riodic boundary conditions in one spatial dimension (cor- 
responding to parallel Casimir plates), then the two- loop 
corrections to the photon polarization tensor lead to a su- 
perluminal propagation of photons. A similarity between 
the effects of Casimir plates and that of finite tempera- 
ture on the propagation of photons has been discussed 
0, @ , and more general conditions that lead to a super- 
luminal propagation of photons have also been identified. 

A similar effect has been found for photons interacting 
with fermions 0. It has been shown that the transverse 
photons coupled to fermions obeying periodic boundary 
conditions are faster than light when propagating per- 
pendicularly and are slower than light when propagat- 
ing parallelly to the compactified dimension. Since our 
fermions obey the usual antiperiodic boundary conditions 
in the compact temporal dimension, one would expect 
no superluminal effects. However, wc demonstrate that 



in a certain region of temperature and chemical poten- 
tial below the symmetry restoration point, massless pi- 
ons propagate faster than light. Moreover, if one of the 
spatial dimensions is compactified, pions will propagate 
superluminally or subluminally, depending on the size of 
the compact dimension and on the boundary conditions. 

We organize the paper as follows. In Scc.[n]we describe 
the model. In Sec. IIIII we calculate the dependence of 
the pion velocity on temperature and chemical potential. 
We present the results and discussion in Sec. IIVI In the 
concluding section, Sec. |Vl we summarize our results. 



II. EFFECTIVE LAGRANGIAN 

Consider the linear sigma model at finite temperature 
and finite baryon density. The thermal bath provides 
a medium which may be assumed to have a homoge- 
neous velocity field. The dynamics of mesons in such 
a medium is described by an effective chirally symmetric 
Lagrangian of the form 

C = -{arf v + bu^u^d^advCT + d^d^) 

-^(a 2 +7 r 2 )-^(a 2 +7r 2 ) 2 
■\-'4>{cvf , d ll + fiu^ + g(cr + iTTrj 5 ))ip , 
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(1) 

where is the velocity of the fluid, and rj^ is the flat 
Lorentzian metric tensor. In the original sigma model 
the fermion field was a nucleon. We consider fermions to 
be constituent quarks [3 El 13 El 13 d 

with addi- 
tional iV c degrees of freedom, "colors", from the SU(iV c ) 
local gauge group of an underlying gauge theory (QCD). 
The parameters a, b, and c depend on the temperature 
T, the chemical potential /x, and the parameters of the 
model to , A, and g, and may be calculated in pertur- 
bation theory. At T = /i = the medium is absent and 
a = c = 1 and & = 0. 

If TOq < 0, chiral symmetry will be spontaneously bro- 
ken. At the classical level, the a field develops a nonvan- 
ishing expectation value (a) = f„. At nonzero tempera- 
ture the expectation value (a) depends on the tempera- 
ture and vanishes at the chiral transition point. Redefin- 
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ing the fields as constant values tt = and a = (a), respectively, we 

obtain the effective Lagrangian in which chiral symmetry 
(a, n) -> (a, tt) + (a'(x),n'(x)), (2) i s spontaneously broken: 

where ir' and a' are quantum fluctuations around the 



-(arf v + bu^u^d^'d^' + -{arf v + bu»u»)d^a 'd v a' 



or „r2 



t _'2 



-(a 12 + Tv' 2 ) 2 + tpici^df, + m F + fiu^ + gW + iriz' 75))^ 



7r —gala + 7r j 



(3) 



r 



At temperatures and chemical potentials below the chiral 
transition line the trilinear coupling and the masses are 
given by 



2Xa 2 , 
0, g' 



m F = ga . 
= Act. 



(4) 



in agreement with the Goldstone theorem. 

The boson kinetic part in Eq. © is split in two terms 
since the pion and sigma self-energies get different finite 
temperature and chemical potential contributions in the 
chiral-symmetry broken phase. Hence, in general, a a 
and b 7^ 6, at temperatures and chemical potentials be- 
low the chiral transition line. The constant c is related 
to the finite temperature and chemical potential contri- 
butions to the fermion self-energy and is irrelevant to the 
calculation of the pion velocity. As we shall shortly see, 
the only quantities relevant to the calculation of the pion 
velocity are the constants a and 6 that enter the pion 
kinetic term in Eq. ©. 

The temperature dependence of the chiral condensate 
a is obtained by minimizing the thermodynamic poten- 
tial VL = —(T/V)\nZ with respect to a at fixed inverse 
temperature (i. At one-loop order, the extremum condi- 
tion reads 



° 2 = ll~^ 
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d 3 p 1 
(2tt) 3 U a 



d 3 p 
(2tt) 3 2^ 



n F (uj F ) 

d 3 p 1 
(2tt) 3 ^ 



where 



us F 



= (P 2 
= (P 2 



n F (u>) 



1 



71 B (uj) 



flu 



(5) 



(0) 



(7) 



The right-hand side of JSJ) depends on a through the 
masses m a and m F given by Solutions to JSJ) are 
implicit functions of T and /i. These equations have been 
derived from the thermodynamic potential in which the 
loop corrections have been neglected [l3|. This approxi- 
mation corresponds to the leading order in the 1/N ex- 
pansion, where N is the number of scalar fields [To ]. In 
our case, N = 4. 



III. PION VELOCITY 

The propagation of pions is governed by the equation 
of motion 

[ (a rf v + b u»u v )] d v ir + V{a, tt, ^)tt = , (9) 

where V is the interaction potential the form of which is 
irrelevant to our consideration. 

In the simplest case of a homogeneous flow, Eq. @ 
reduces to the wave equation 

(d 2 - v 2 A + —V)n = 0, (10) 
a 

where the quantity v is the pion velocity given by 

" 2 =H)"- <") 

As we shall shortly demonstrate, the constants a and b 
may be derived from the finite-temperature perturbation 
expansion of the pion self-energy. 

The pion velocity in a sigma model at finite tempera- 
ture has been calculated at one-loop level by Pisarski and 
Tytgat in the low-temperature approximation [2j and by 
Son and Stephanov for temperatures close to the chiral 
transition point 0, 13 ■ It has been found that the pion 
velocity vanishes as one approaches the critical temper- 
ature. Here we analyze the whole range of temperatures 
in the chiral symmetry broken phase. 

Consider the pion self-energy T) in the limit when 
the external momentum q approaches 0. The renormal- 
ized inverse pion propagator may, in the limit q — > 0, be 
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expressed in the form 

+ ..., ' ' (12) 

where the ellipsis denotes the terms of higher order in 
q^ . The q^ independent term of the self-energy absorbs 
in the renormalized pion mass, which is equal to zero in 
the chiral symmetry broken phase. The subtracted T=0 
term has been absorbed in the wave function rcnormal- 
ization factor Z n . By comparing this equation with the 
inverse pion propagator derived directly from the effec- 
tive Lagrangian (pj) 



A" 



(a + b)ql 



aq 2 , 



(13) 



we can express the parameters a and 6, and hence the 
pion velocity, in terms of the second derivatives of E(g, T) 
evaluated at q^ — 0. At one-loop level the diagrams 
that give a nontrivial q-dependence of E are the bubble 
diagrams. Subtracting the T=0 term one finds 



E(g) = E(g, T) - E(g, 0) = E B + E F 



(14) 



with the contribution of the bosonic and fermionic loops 
given by 



Es(g) = -V 2 / -/^JL _ L <{ [ nB (u3^) + n B {u> a ,q)] 



(2tt) 3 2w w 2w crjg 
+ [n B (uir) - n B (uj a ^)} 



1 



1 



qo 



qo 



(15) 



E F (q) = -8N c g 



d 3 p 



-qo^F.q + q{p + q) 



(2ir) 3 2ujF2ujF,q I e /3 (" F .9 *0 + 1 \ojf + &F,q — qo up — up, q + qo 
q 2 - q uj F + q(p + q) ( 1 1 



+ 



e /3(^F M) + 1 \WF,g + WF + 90 tOp_ q - UJp - q 

qo^F,q + g(p + g) / i ^ i 

e /3(^F,,+M) ^CJF + U>F,q + qo tOF - ^F,q ~ qo 

-q 2 + qauJF + q{p + q) ( I i 



e [3(uj F +n) _|_ I \u} F ,q+UF-qo UF,q — WF+qo 



(16) 



where u a>q = [(p — q) 2 + m 2 } 1 / 2 , ujp.q = [(p + q) 2 + mp] 1 / 2 . A straightforward evaluation of the second derivatives 
of E(g) at q^ = yields 



a = 1 + a B + ap ■ 



(17) 



6 = 6, 



(18) 



with 



a B 



d 3 p 
(2Vp 



4w T 



n B (uJa) 1 u)\ ( n B (u> v ) _ 71_b(u ct ) 

4Wcr 3mJ \ CJtt U) a 



(19) 



16/ 



d 3 p 
(2^)3 



w v n B (uj^) u a n B (uj a ) 1 w 2 / 'n B (u n ) n B (u a ) 



r 



3mJ \ ^ 



(20) 



a F = N c g 2 



d 3 p tif(u>f) 
(2ir) 3 p 2 u)F 



(21) 



b F = -N c g 2 



d 3 p m 2 F uf^f) 



(22) 



(2tt) 3 p 2 lu 3 f ■ 
The subscripts B and F denote the contributions of the 
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boson and fermion loops, respectively. The sign of the 
fcrmion contribution in the last equation is crucial. The 
pion velocity v given by Eq. 1)11(1 will become larger than 
unity when b < 0, i.e., when the negative fermionic part 
exceeds the positive bosonic part. 

The results JE) and JTHJ) with fT5 |l -ip2 )l can easily 
be extended to include the sigma model at T = and 
fi = in 3+1 dimensional space-time with one spatial 
dimension, e.g., in the z-direction, compactified to the 
size L = (3. In this case, the effective Lagrangian may 
be written in the form of Eq. Q with b replaced by —6 
and the velocity m m replaced by a spacelike vector 
parallel to the compactified dimension and normalized 
as n^n^ = — 1. The inverse propagator is now given by 



A" 



a (<7o 



ql)-{a + b)ql 



(23) 



and hence the velocity in the compact direction precisely 
equals the inverse of v given by Hill) . Clearly, the veloc- 
ity in the subspace orthogonal to the compact dimension 
is equal to the velocity of light. Recall that the exten- 
sion of field theory to nonzero temperature is equivalent 
to a compactification of the time coordinate. Now, the 
compactified spatial coordinate and the time coordinate 
exchange their roles. However, since there is no time- 
ordering restriction in spatial coordinates, fermions may 
be chosen to be periodic or antiperiodic in the compact 
dimension. Hence, in the case of antiperiodic fermions, 
all the calculations of the self-energy are exactly the same 
as in the finite temperature case with T = l/L. The case 
of periodic fermions is similar, but one has to take into 
account that the replacement of the antiperiodic bound- 
ary condition at \x — by the periodic one modifies the 
/3-dependent part of the loop sum over 4-momenta so that 
the Fermi-Dirac function 1/(1 + e' 3 ") in uf is replaced 
by 1/(1 - e^). In this way the pion velocity parallel to 
the compact direction may be written as 



«jj = l+^, (24) 
where a' — a and b' = b for antiperiodic and 

a = 1 + clb — ap , (25) 



b' = b B - b 
for periodic fermions, with 

_> , d 3 p n B (uj F ) 



(26) 



a F = 2N c g' 



b F = -2N c g 2 



(2ir) 3 p 2 u>F ' 



d 3 p m F n B {uJF) 
(2tt) 3 p 2 uj 3 f 



(27) 



(28) 



Note that in the case of periodic fermions Eq. (J5J 
also modifies, because np(bjp) must be replaced by 
— 2n_e(w_F). Depending on the choice of the coupling con- 
stants g and A, the symmetry breaking and restoration 
pattern may be quite different. In particular, the temper- 
ature T' c at which sigma vanishes and which approximates 
the critical temperature |l3j |. is now given by 



GA 



-2g 2 N c + 3A 



(29) 



Clearly, for A < (2/3)g 2 N c there is no restoration of chi- 
ral symmetry. 



IV. RESULTS AND DISCUSSION 

In order to proceed with a quantitative analysis, we 
have to choose the input parameters from particle physics 
phenomenology. For the constituent quark mass we take 
mp = 340 MeV. The coupling g is then fixed by Q 
g = viif I ' f-K- The Particle Data Group gives a rather wide 
range 400-1200 MeV 01 of the sigma meson masses. We 
shall shortly see that the analysis of the T = Q, /j, ^ 
case further restricts this available range. 

At T = the extremum condition JSJ) for nonnegative 
fi and N c = 3 reads 



\x\J n 2 — g 2 a 2 + g 2 a 2 In ■ 



.90" 



+ v V 2 - g 2 & 2 



which determines the extremum of the thermodynamic potential 



Ofn- ,A A - 4 A f 2 - 2 tit Z gg ) 



/Va* 2 -9 2 o 2 {2^ 2 - bg 2 a 2 ) + 3<?V 4 In £ 
I 



+ yV 2 - 9 2 cr 2 

go- 



(30) 



(31) 



The chemical potential /i c at which chiral symmetry gets restored is given by the condition $1(0, fi c ) = £l(a,fi c ). 
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FIG. 1: Pion velocity squared as a function of temperature and chemical potential. The dashed line corresponds to a 
thermodynamically unstable solution. T and fi are in MeV. 



For \i c < gf n , the solution of is a = f^, so lp?0|) 
leads to 



/'<-: 



(32) 



Since // c is the threshold where extended nuclear matter 
forms, a reasonable assumption, confirmed by the strong- 
coupling QCD analysis is 



(33) 



It may be easily verified that with this assumption the 
pion velocity is constant and equal to light velocity up to 
the transition point /i c , where it drops to zero. 

It is remarkable that Eq. I|33|l yields an upper bound 
on the sigma mass 



2m% 



(34) 



If we take / w = 92.4 MeV, m F = 340 MeV, and saturate 
the bound, we find m a — 796.5 MeV. In the following 
numerical analysis we take this value as an input param- 
eter. Any other choice below this value would not alter 
our qualitative picture. 



Next we analyze the case of nonzero temperature. In 
Fig.^we plot the pion velocity v as a function of temper- 
ature for various fixed fi (upper plots) and as a function of 
chemical potential for various fixed T (lower plots). For 
each T and ji the chiral condensate a is calculated by 
solving Eq. JSJ numerically. The dashed line represents 
the velocity obtained with thermodynamically unstable 
solutions for a. This line follows a similar backwards 
going line in the plot for a 13], which is typical of the 
first order phase transition. The actual phase transition 
takes place at the point T c (for a given [i = fj, c ) or at 
\x c (for a given T = T c ) where the two minima of the 
thermodynamic potential at a = and <j(T c ,(i c ) are lev- 
eled 01 . At the critical point, the pion velocity drops 
abruptly to zero. We note that there is always a region 
of temperatures or chemical potentials below the criti- 
cal point where the pion velocity becomes superluminal. 
The temperature T = 50 MeV and the chemical potential 
fi = 160 MeV are chosen to represent a typical behavior 
below the critical temperature and the critical chemical 
potential, respectively. The value T = 3 MeV is cho- 
sen to show the low-temperature behavior near the point 
T = which we have discussed above. If we increase 



6 



the temperature, starting from the point (T = 0, fj, = 0) 
where v = 1, the velocity at the beginning drops slightly 
below 1 and at some temperature Tsl of the order of 50 
MeV, it becomes superluminal. 

It is important to check how sensitive the onset of su- 
perluminal propagation is to the variation of the param- 
eters of the model. In Fig. we plot the superluminal 
onset temperature TgL as a function of the fermion-boson 
coupling g keeping the boson self-coupling A fixed. Con- 
trary to what one would naively expect, the superluminal 
onset shifts to larger temperatures with increasing g. The 
reason is that in this temperature range the contribution 
of fermions (|22|l . relative to that of bosons (|20ll . decreases 
with increasing g owing to the fermion mass dependence 
of the distribution function. 

The behavior of the pion velocity near the critical tem- 
perature should be analyzed with special care. It is well 
known that the phase transition in the SU(2)xSU(2) lin- 
ear sigma model at fi = is of second order [2Cj and 
that the one-loop approximation violates this prediction 
by producing a weak first-order transition l2lj | . Cal- 
culations of the chiral condensate in the chiral perturba- 
tion theory j2^ show the deviation from nonperturbative 
renormalization group calculations |23j at temperatures 
of about 0.5 T c . As our results depicted in Fig. ^ (upper 
left plot) show that the superluminal pion propagation 
starts at temperatures of about 50 MeV, one would per- 
haps be inclined to think that the superluminal effect 
was an artifact of the one-loop approximation. However, 
as we have already pointed out, the superluminal pion 
velocity is a consequence of the negative fermion loop 
contribution. The proportion of the fermion part to that 
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FIG. 2: Superluminal onset temperature as a function of the 
fermion-boson coupling g. 



2 






i 
/ 








/ 








/ 








/ 








/ 








/ 


- 






/ 








/ 








/ 








i / 








v 

i 

i 
i 






1/L [MeV] 



50 100 150 



FIG. 3: Pion velocity squared as a function of inverse com- 
pactification length with antiperiodic (solid line) and periodic 
(dot-dashed line) fermions. The dashed line corresponds to a 
physically unstable solution. 



of the boson depends on the parameters of the model. 
For example, a decreasing of g with other free param- 
eters fixed, increases the negative fermion contribution 
given by Eq. (|22ll . which in turn yields a superluminal 
propagation even at temperatures much below 50 MeV 
(Fig. |2Jl. Hence, we believe that the onset of superluminal 
propagation coincides only accidentally with the window 
where perturbation theory starts to deviate from non- 
perturbative calculations. For a firm judgment whether 
a superluminal pion propagation is a genuine effect or 
an artifact of the perturbation expansion one should go 
beyond the one-loop approximation or apply a nonper- 
turbative technique such as Monte Carlo calculations on 
the lattice or the 1/N expansion of the 0(N + l)/0(N) 
nonlinear sigma model with fermions. 

Next we discuss the case T = and /i = with one 
spatial dimension compactified to the size L = [3. The 
pion velocity v\\ parallel to the compact direction is plot- 
ted in Fig.[3]for both periodic and antiperiodic fermions. 
In the case of antiperiodic fermions, the compactification 
size plays the role of inverse temperature and chiral sym- 
metry gets restored at the critical size L c = l/T c . In this 
case, the velocity plot is just the inverse of the first plot 
in Fig. □ 

For periodic fermions we find the minimum of the ther- 
modynamic potential corresponding to a nonzero sigma 
to be always below the minimum corresponding to a = 0. 
In this case, the right-hand side of l|29|l is negative and 
no restoration of chiral symmetry takes place. Owing to 
the opposite sign of the fermion contribution the veloc- 
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ity uii is always superluminal, monotonously approaching 
the velocity of light as the size of the compactification L 
approaches infinity. 

To gain a physical insight into why the pions become 
superluminal it is worthwhile making a comparison with 
the similar effects found for photons propagating in a 
medium and in a Casimir vacuum. In QED at either fi- 
nite temperature or finite boundaries the genuine change 
of the speed of massless photons is a two loop effect re- 
lated to the coherent light-by-light scattering 0,Hi- How- 
ever, a one loop superluminal effect has been found in 
QED with one space dimension compactified |9j . Modifi- 
cations of the vacuum that change the population of real 
and virtual particles introduce coherent scattering which 
decreases or increases the speed of massless photons. In 
a Casimir vacuum some of the virtual modes are elim- 
inated and consequently their would be scattering 
As a result the speed of massless photons is increased. 
Roughly speaking, the negative Casimir energy density 
results in massless photons being superluminal. In our 
case, a superluminal propagation is, formally, a conse- 
quence of the negative fermion one-loop contribution to 
the q 2 coefficient in the inverse pion propagator It 
is conceivable that, for the parameters close to the phase 
transition point, the specific interaction between the chi- 
ral fields and fermions reduces the effective vacuum en- 
ergy density and as a result the massless pions become 
superluminal. 

V. SUMMARY AND CONCLUSION 

We have analyzed the propagation of massless pions at 
nonzero temperature and nonzero finite baryon density 
in a sigma model with two quark flavors. By calculating 
the pion dispersion relation at one-loop order we have 
shown that pions in the presence of fermions become su- 
perluminal in a certain range of temperatures and baryon 
chemical potential. 

Furthermore, we have studied the case when one of the 
spatial dimensions is compactified with fermions obeying 
periodic or antiperiodic boundary conditions. Restrict- 
ing attention to T = p, = 0, we have calculated the pion 
velocity v\\ along the compact direction. We have found 



that for antiperiodic fermions, pions will propagate su- 
pcrluminally or subluminally, depending on the size of 
the compact dimension. With periodic fermions, the ve- 
locity v\ I is always larger than the velocity of light. 

A superluminal propagation of massless particles may 
naively seem to contradict special relativity and Lorentz 
invariance. The most disturbing consequence of a su- 
perluminal propagation would be an apparent violation 
of causality. However, it has been convincingly argued 
that superluminal effects realized in quantum field theo- 
ries with a nontrivial vacuum do not lead to causal para- 
doxes [24J. The basic argument goes as follows: once 
the conditions that determine the vacuum fluctuations 
are fixed, the propagation velocity in a given reference 
frame is unique. This implies that it is not possible to 
send signals both forwards and backwards in time within 
one reference frame. In other words, the causal loops are 
forbidden. 

It is important to bear in mind that we have con- 
sidered an idealized situation when there is no explicit 
chiral symmetry breaking of the original Lagrangian l[TJl 
and pions are exactly massless in the broken symmetry 
phase. In reality, chiral symmetry is explicitly broken 
owing to nonvanishing current quark masses. To make 
the pions massive we could have, as usual, introduced 
an explicit chiral symmetry breaking term in the origi- 
nal Lagrangian. In that case the quantity v, defined in 
(fTTf) and calculated in the limit q — > 0, would no longer 
have the meaning of the pion velocity. Instead, given the 
dispersion relation u> = \J v' 2 q 2 + to 2 , with both v and 
to being T and \i dependent, one would define a group 
velocity v g — dcu/dq. It may easily be verified that for 
q <c to the quantity v g is below one even if v > 1. Of 
course, the situation could change for q ~ to. That would 
involve a calculation of v at non-zero momenta which is 
beyond the scope of the present paper. 
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